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1.  INTRODUCTION 


In  the  preceding  paper  (Jensen,  Adams,  and  Chabalowski,  to  be  published)  (Referred  to  as  Report 
I)  we  reported  ab  initio  calculations  predicting  the  magnetic  circular  dichroism  (MCD)  spectra  for  the 
C< — X  transition  in  acetylene.  In  this  second  paper  on  MCD,  we  report  on  a  new  application  of  a 
recently  implemented  method  for  studying  electronic  spectral  properties  which  exist  due  to  interactions 
(perturbations)  of  the  electronic  state  TA  with  a  manifold  of  electronic  states  (T/s,  with  J*A).  This 
type  of  intensity  borrowing  is  fundamental  to  the  ®0-term  found  in  the  general  expression  for  the  MCD 
intensity.  In  this  study,  we  apply  this  new  methodology  to  the  a  priori  calculation  of  MCD  spectra, 
again  in  the  acetylene  molecule,  but  for  the  spectral  region  185-170  nm.  The  method  is  based  on  first- 
order  perturbation  theory  and  has  been  used  successfully  to  study  spin-orbit  and  spin-spin  interactions 
in  other  molecules  (Havriliak  and  Yarkony  1985;  Chabalowski  et  al.  1989). 

Gedanken  and  Schnepp  (1976)  have  published  the  experimental  electronic  absorption  and  MCD 
spectra  for  acetylene  in  the  spectral  region  190-170  nm.  They  did  not  attempt  to  make  an  assignment 
or  further  explain  the  structure  of  a  band  system  occurring  in  this  region.  Other  experimental  studies 
(Ingold  and  King  1953;  King  and  Ingold  1952;  Lassettre  et  al.  1968;  Dance  and  Walker  1973;  Foo  and 
Innes  1973)  have  added  strong  evidence  to  the  onset  of  a  new  band  system  near  185  nm,  with  the 
primary  assignment  of  these  spectra  to  the  singlet  Bt-X  transition,  with  B  taken  as  the  l'Bu  state  in  the 
Qh  point  group  for  a  trans-bent  (planar)  conformation.  Figure  1  shows  the  symmetry  of  the  states 
arising  from  splitting  the  doubly  degenerate  l'A^  state  (linear  conformation)  by  bending  into  the  trans 
and  cis-conformations.  The  spectrum  beginning  at  185  nm  is  reported  to  be  "complicated"  and 
"diffuse"  (Gedanken  and  Schnepp  1976),  due  perhaps  to  some  overlap  by  other  electronic  transitions 
as  well  as  vibrational  coupling  of  the  B  to  other  electronic  states  (Peric,  Buenker,  and  Peyerimhoff 
1984).  This  is  consistent  with  the  experiment  of  Foo  and  Innes  (1973)  who  have  assigned  high-lying 
vibrational  bands  of  the  A(lAu)<-X  transition  at  energies  near  the  0-0  of  the  B<— X 

Theoretical  studies  have  been  canned  out  on  several  singlet  electronic  states  in  acetylene  (Peric, 
Peyerimhoff,  and  Buenker  1985,  1987;  Peric,  Buenker,  and  Peyerimhoff  1984;  Kammer  1970,  1974; 
Lischka  and  Karpfen  1986;  So,  Wetmore,  and  Schaefer  1980)  and  segments  of  the  theoretically 
predicted  potential  energy  surfaces  (PES)  were  presented  for  these  states  (Peric,  Peyerimhoff,  and 
Buenker  1987;  Kammer  1970).  The  results  of  the  theoretical  work  support  the  assignment  of  at  least 
part  of  the  intensity  in  the  spectral  region  185-170  nm  to  the  l'Bu  state,  and  this  assignment  now 
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seems  well  accepted.  The  theory  also  suggests  that  the  llB2«-X  transition  could  fall  within  this 
spectral  region,  causing  part  of  the  irregularities  seen  in  the  spectra. 

Since  only  singlet  electronic  states  are  treated  in  this  study,  the  spin  multiplicity  will  be  dropped 
from  the  state  labels.  The  B  (Bu)  state  is  the  partner  to  the  2AU  state  which  arises  from  splitting  the 
lA^  state  by  a  planar  trans-bending  motion  of  the  hydrogens  (see  Figure  1).  The  B  state  shows  a 
minimum  in  its  PES  near  a  140°  bend  angle  of  the  C-C-H  with  the  linear  conformation  taken  as  180°. 
Because  the  B  state  is  not  degenerate  near  its  equilibrium  point,  it  will  be  assumed  that  the  MCD 
intensity  arises  from  the  ®0-term  which  requires  the  use  of  non-degenerate  perturbation  theory  to 
account  for  the  intensity  borrowing  from  other  electronic  states  through  magnetic  dipole  interactions. 
The  B«— X  transition  will,  therefore,  be  treated  using  the  newly  implemented  first-order  perturbation 
techniques  mentioned  above.  The  2AU  state  is  predicted  to  be  linear  (Peric,  Peycrimhoff,  and  Buenker 
1987;  Peric,  Buenker,  and  Peyerimhoff  1984;  Rammer  1970;  So,  Wetmore,  and  Schaefer  1980)  if  the 
molecule  is  kept  planar.  However,  earlier  configuration  interaction  (Cl)  calculations  show  that  the  2\ 
state  correlates  with  the  3'A  state  in  Q  symmetry,  and  the  3'A  has  a  bent  and  twisted  (about  the  C-C 
bond)  minimum  energy  structure  (Peric,  Peyerimhoff,  and  Buenker  1987).  So  excitation  into  the  2A^ 
state  should  result  in  a  bent  and  non-planar  structure  of  Q  symmetry. 

Transition  energies,  absorption,  and  MCD  intensities  will  also  be  calculated  for  the  1B2<— Xband 
system  in  order  to  consider  the  possibility  of  this  transition  contributing  to  the  spectra  in  the 
185-170  nm  region.  The  1B2  arises  from  the  cis-bending  which  splits  the  l/\,  (Figure  1)  into  1B2 
and  2Aj  in  the  Q,  point  group,  with  the  1B2  predicted  to  have  a  non-linear  (but  planar)  equilibrium 
geometry  at  a  bend  angle  near  130*"  (Peric,  Peyerimhoff,  and  Buenker  1987;  Peric,  Buenker,  and 
Peyerimhoff  1984;  Rammer  1970;  So,  Wetmore,  and  Schaefer  1980).  The  2A^  (like  the  2AU  of  the 
trans-bent)  correlates  with  the  3‘A  state  in  Q  and  the  same  comments  apply  to  its  minimum  energy 
structure  (vide  supra).  For  a  comprehensive  review  of  the  singlet  electronic  states  in  acetylene,  the 
reader  is  referred  to  the  paper  by  Peric,  Peyerimhoff,  and  Buenker  (1987). 

The  ab  initio  wavcfunctions  will  be  developed  from  large  atomic  Gaussian  basis  sets  and  state- 
averaged  multiconfiguration  self-consistent  field  (MCSCF)/CI. 
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2.  METHODS 


State-averaged,  complete  active  space  MCSCF  (SA-CASSCF)  was  performed  for  several  electronic 
states  at  each  point  on  the  PESs  and  the  resulting  molecular  orbitals  (MO)  were  used  as  expansion 
vectors  in  the  CIs.  The  Cl  wavefunctions  were  constructed  in  the  space  of  all  single  and  double 
excitations  from  a  set  of  reference  configuration  state  functions  (CSF).  Information  on  the  SA- 
CASSCF  and  Cl  procedures  used  to  obtain  the  zeroth-order  wavefunctions  was  presented  in  Report  I 
and  will  not  be  repeated  here.  Details  of  their  applications  to  this  study  will  be  presented  in  the 
appropriate  sections.  For  the  sake  of  clarity,  the  general  MCD  equations  and  the  details  of  the  MCD 
^,-term  are  presented  here,  as  are  the  perturbation  theory  treatment  of  the  3,-term,  and  the  equations 
needed  to  fit  the  MCD  band  maxima,  denoted  Ae““. 

2.1  MCD  Equations.  The  MCD  intensity  is  given  by  the  difference  between  the  decadic  molar 
extinction  coefficients  for  left  and  right  circularly  polarized  light,  i.e.,  Ae  =  tL  -tK,  (Stephens  1976; 

Piepho  and  Schatz  1983)  and  for  a  transition  involving  a  non-degenerate  ground  state,  X,  and  a 
possible  degenerate  excited  state,  I,  we  get 


££  =  653.2 
u 


^  A'  (^)  +  B°/]‘// 


(1) 


where  i>  is  energy,  f  a  line  shape  function,  H  the  magnetic  field  in  Gauss,  and  pB  the  Bohr  magneton 
(.4669xlO_4cm'1/Gauss).  and  Q  are  given  by 


A  -  x  E  [<4 1? I/,,* •(<*!«  |/,>  x  </j,|m|x>)l 

J  X,X’ 


(2) 
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E 
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%  ~EX 


__  ^  — 

+  £  *  (<X|m  |/i>x</4  |m  |X>) 

Jk*h  EJt'\ 


(3) 


with  ji  and  in  being  the  magnetic  and  electric  dipole  moment  operators.  Equations  2  and  3  take  into 
account  the  fact  that  the  molecules  are  randomly  oriented.  The  units  are  cgs  with  energies  in  cm'1, 
electric  dipole  moments  in  Debye,  and  the  magnetic  field  in  Gauss.  The  operator,  m,  is  chosen  in  the 
length  form,  and  the  magnetic  dipole  moment  operator,  p*.  is  related  to  the  orbital  angular  momentum 
operator,  if,  by  ]x=(iP.)L.  The  two  operators  are  identical  (in  atomic  units)  except  for  the  factor  of 
1/2* 


As  Equation  2  is  written,  one  sees  that  the  !\  term  is  derived  for  a  truly  degenerate  excited  state, 

I ,  and  that  a  MCD  band  attributed  to  such  a  state  will  exhibit  a  characteristic  derivative  line  shape.  In 


*  The  actual  expression  representing  the  Ae/E  calculated  in  this  study  is  given  by 


t-'o,525KI)‘sVJ-" 

which  differs  from  Equation  1  in  the  units  conversion  due  to  the  substitution  of  the  orbital  angular  momentum  operator  for 
the  magnetic  dipole  moment  operator  and  carrying  through  the  multiplication  by  p,,  the  Bohr  magneton.  The  definitions  of 
A,  and  Ot.  then  become 


Al  =  ^  E  ,[<rx\L\rx,>-(<i\m \I\>  x«/4,|w|*>)] 

3  A,A 


E 

Jk*X 


<yt|Z|x> 

Ejk  ~Ex 


•  (  <  X  |  m  |  ik  >  x  <  |  m  |  JK  >) 


♦  £  -iliii-  -(<Jf|m|/;>  X  </Jm|X>) 

Ejk  ~  £/i 


where  use  has  been  made  of  the  relationship  p=(l/2)L  .  The  matrix  elements, 

<IX  \L  |/x,  >,  <Jt\L  |X>,  and  </x  \L  |yt>  come  directly  from  the  quantum  chemical  calculations. 
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contrast,  the  ^,-terai  can,  in  general,  be  non-zero  for  any  electronic  transition  and  will  have  the  simple 
line  shape  of  an  absoiption  band.  These  line  shape  differences  are  often  used  to  determine  if  a  band 
involves  a  transition  to  a  degenerate  excited  state.  It  should  be  pointed  out  that  the  %term  clearly 
contains  the  form  of  a  first-order  perturbation;  i.e.,  X  and  I  states  are  perturbed  by  a  manifold  of 
J  states  via  the  magnetic  dipole  moment  operator.  For  perturbations  to  the  I  state: 


E 


<4lnl4> 


•</„ 


(4) 


The  reader  is  referred  to  Stephens  (1976)  and  Piepho  and  Schatz  (1983)  for  detailed  discussions  of 
MCD  theory,  applications,  and  band  shapes. 

In  order  to  relate  the  observed  spectra  to  our  theoretical  results,  we  use  a  normalized  Gaussian  as 
our  line  shape  function: 


/=  _L_e-('>-'5,)2/A2  (5) 

A  (IT)1/2 


^  -2(v-5,)  e_ (v - u0)2/A2  .  (6) 

A3  (II)W 


The  A  is  one-half  the  bandwidth  measured  at  the  band  height  of  (fmc/e)  and  generally  taken  from  the 
absorption  spectrum.  The  transition  energy  in  wavenumbers  is  v  and  \>0  is  the  transition  energy  at  the 
band  maximum.  Due  to  the  lack  of  a  baseline  in  the  absorption  spectrum  of  Gedanken  and  Schnepp, 
rough  estimates  of  the  bandwidths  were  taken  from  the  MCD  spectrum  of  Gedanken  and  Schnepp. 


6 


(7) 


The  oscillator  strength  for  absorption  from  state  X  to  state  I  is  defined  as 

/*,«  |A£(au)|<X|mV>|2 

where  AE  is  the  transition  energy  and  all  quantities  are  in  atomic  units. 

2.2  Perturbed  Wavefunctions.  The  perturbed  wavefunction  for  state  %  is  given  through  first- 
order  in  perturbation  theory  by 

*¥,  «  +  'P1/  (g) 


with  *Fj  being  the  zeroth-order  approximation.  The  usual  spectral  representation  for  the  first-order 
correction.  *F},  due  to  magnetic  dipole  interactions  would  be 

e  (9) 

ui  (e,°'e;) 

The  summation  over  the  L  electronic  states  is,  in  principle,  infinite.  One  often  used  approach  to 
solving  for  is  to  calculate  explicitly  the  wavefunctions  for  a  relatively  small  number  of  excited 
states,  thereby  drastically  truncating  L.  This  might  cause  one  to  miss  important  contributions  to  4'1, 
from  the  omitted  states. 


Within  a  given  CSF  space,  this  "omitted  states"  problem  is  eliminated  by  solving  for  4*',  directly 
from  (Haviiliak  and  Yaikony  1985) 

-p'F®.  00) 

Equation  10  can  be  transformed  into  matrix  form  as 

{.H°-E)V' •-w:' 

(11) 
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where  it  must  be  emphasized  that  Hand  (i  matrices  with  elements  formed  over  CSFs,  not  over 
eigenstates.  The  vectors  V  and^?  are  defined  as  the  coefficients  for  the  first-  and  zenoth-order  parts  of 


=  EC/fiW  (12a) 

i 

-  EV/V/k*).  (12b) 

s 

The  k  and  k1  label  the  spatial  symmetries  to  which  the  CSFs  belong,  and,  in  general,  k*k\ 

Equation  1 1  represents  a  large  set  of  linear  inhomogeneous  equations  which  are  solved  to  obtain  V  by 
a  variant  of  the  method  suggested  by  Pople,  et  al.  (1979). 

3.  DETAILS  OF  CALCULATIONS 

The  basis  set  used  is  identical  to  that  described  in  Report  I,  and  consists  of  the  primitives  (1  Is,  7p) 
contracted  to  (5s,  4p)  on  the  carbons  and  (6s)  contracted  to  (3s)  on  the  hydrogens.  The  carbon  atoms 
were  augmented  with  two  uncontracted  d-type  polarization  functions  and  the  hydrogens  each  received 
one  uncontracted  p-type  function.  To  the  center  of  the  C-C  bond  was  added  a  set  of  (2s,  2p,  Id) 
uncontracted  diffuse  functions,  thus  giving  a  total  of  84  atomic  orbitals  (AO),  with  the  details  given  in 
Table  1  of  Report  I.  The  details  of  the  SA-CASSCF  and  Cl  calculations  will  be  presented  in  the 
sections  discussing  the  electronic  states. 

Peric,  Peyerimhoff,  and  Buenker  (1987)  have  discussed  many  of  the  difficulties  associated  with 
studying  the  vibrational  motions  of  the  excited  electronic  states  in  acetylene,  and  particularly  the  B 
state.  Not  the  least  of  these  difficulties  is  the  probability  of  strong  couplings  amongst  the  bending  and 
twisting  modes.  However,  a  treatment  of  the  nuclear  motion  is  required  here  because  of  a  large 
difference  in  the  equilibrium  structures  of  the  ground  state  and  the  B  and  IB,  states.  A  semi- 
quantitauve  treatment  of  the  vibronic  problem  included  in  this  study  should  be  sufficient  to  address  the 
electronic  origin  of  the  spectral  features  in  the  185-170  ran  range  which  is  the  goal  of  this  work. 
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In  this  work,  we  employ  a  simplified  treatment  of  the  nuclear  motion,  ignoring  all  vibrational 
mode  couplings  (excepting  the  2au-B  over  p*.  to  be  discussed  later)  and  keeping  the  C-C  and  C-H  bond 
lengths  fixed  at  the  ground  state  values  of  R(C-C)=1.208A  and  R(C-H)=1.058A20  for  the  calculation  of 
transition  moments.  This  choice  of  bond  lengths  was  based  upon  the  fact  that  we  are  predicting  an 
absorption  intensity  from  the  ground  state.  The  PESs  for  two  states  were  then  recalculated  using  the 
bond  lengths  determined  by  Peric,  Peyerimhoff,  and  Bucnker  for  the  minimum  energy  structure  of 
each  state,  these  being  the  B[R(C-C)=1.349A,  R(C-H)=1.07lA]  and  lB2[R(C-C)=1.35lA, 
R,(C-H)=1.104A].  These  PESs  were  then  used  in  calculating  a  more  appropriate  set  of  vibrational 
wavefunctions.  The  PES  for  the  2Att  was  uniformly  lowered  in  energy  by  the  energy  difference  for 
the  B  state  calculated  at  180°  using  these  two  different  sets  of  bond  lengths.  This  was  necessary  to 
facilitate  a  more  accurate  coupling  of  the  B  and  2AU  states. 

The  bending  motions  were  described  as  a  harmonic  variation  in  the  bending  coordinate,  0,  with  0 
given  in  radians.  Due  to  the  symmetry  of  the  bending  modes,  the  moments  between  the  V"=0  (gerade) 
and  all  the  odd  numbered  v'  (ungerade)  levels  must  integrate  to  zero  over  the  electric  dipole  moment 
matrix  elements.  So  the  calculated  vibrational  levels  describing  the  zero  intensity  moments  will  be 
removed  and  the  remaining  even  states  relabeled  consecutively  to  correspond  with  the  experimental 
notatioa 

Within  the  context  of  a  single  mode  approximation  it  is  necessary  to  assign  a  reduced  mass  for  the 
mode  in  question.  The  reduced  mass  was  obtained  by  fitting  the  calculated  fundamental  bending 
frequencies  vt  (trans)  and  us  (cis)  to  experiment  for  the  X  and  B,  or  to  the  calculated  harmonic  values 
of  Lischka  and  Karpfen  (1986)  for  the  2A„  state.  This  treatment  is  expected  to  be  more  reliable  for 
the  lower  vibrational  levels.  The  reduced  masses  for  the  trans-bending  mode  are  calculated  to  be: 

(X  :1.0),  (B  :0.73),  and  (2A«  :0.57).  The  reduced  masses  for  the  cis-bending  mode  are:  (X  :1.8), 

(1B2  : 1 .3),  and  (2A2:0.635). 

Henceforth,  all  quantities  shall  be  reported  in  atomic  units  unless  otherwise  stated.  The  trans-bent 
and  cis-bent  conformations  are  treated  in  the  point  groups  C^,  and  Q,,,  respectively,  and  the 
correspondence  of  the  irreducible  representations  (IRREP)  amongst  the  three  symmetries,  D^,  C*,,  and 
Q,,  is  given  in  Table  1.  corresponds  to  the  linear  molecule  and  is  included  for  explanatory 
purposes.  The  molecule  lies  in  the  xy  plane  in  Cy,  with  the  C-C  bond  collinear  to  the  x-axis,  and  in 
Cjv  the  molecule  is  in  the  yz  plane  with  the  C-C  collinear  to  the  y-axis. 
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Table  1.  Correspondence  of  Irreducible  Representations  for  Acetylene  in  the  Trans-bent  C^,  Linear 
D^,  and  Cis-bent  Qv  Conformations 


Qh 

Dooh 

Cjy 

A, 

A, 

A„ 

A2 

Au,  Bu 

Au,  Bu 

A  ,  B 

A 

M 

nu 

nf 

a2,  b2 

A],  Bj 

a2,  b2 

r  g 

8 

A,,  B, 

A..  Bg 

A, 

3.1  The  R— X  Transition:  Trans-bend. 

3.1.1  Zeroth-order  State  Descriptions.  The  trans-bend  B  (B„)  state,  which  arises  from  the 
degenerate  1A„,  is  represented  by  the  CSFs  ( lo* \c\la\2a\ )  30*17^  Itt^  Itt^,  and 

(.  .  .)  3o*  lrc^  1  ni,  1  jc .  Its  degenerate  partner  is  the  linear  2A„  and  is  given  by  the  CSFs 
(.  .  .)  30^1^  ljj,l  jtj,  and  (.  .  .)  30,171^171^17^.  The  active  orbitals  used  in  the  SA-CASSCF  are 
the  3<r<(4a<),  171^(5^),  ljt>Q((36ii)tljtiii(/aii),l7C4J(16p  with  the  4rcu  electrons  distributed  amongst  these 
MOs.  The  3o,(4 ag)  MO  is  a  Rydberg  orbital.  Six  electronic  states  were  then  averaged  including 
1  Ag, \2Bu,\,2,3Am.  The  weighting  scheme  used  in  the  averaging  is  wk=(l, 1,1,1, 1,1),  respectively. 
The  SA-CASSCF  was  calculated  at  each  point  along  the  bending  PES,  which  ranged  from  180° 
(linear)  to  110°  at  10°  increments,  and  another  point  at  90°. 

To  test  the  effects  of  the  choice  of  weights  on  the  transition  moments,  a  test  was  run  at  a  bend 
angle  of  150°  using  the  new  weights  wk=(2,l,l,l,l,l).  Only  two  of  the  transition  moments  changed 
more  than  4%,  and  these  two  moments  are  relatively  small  and  contribute  to  the  total  MCD  intensity 
in  such  a  way  as  to  have  changed  the  final  intensity  by  less  than  2%.  These  size  fluctuations  in  the 
moments  are  quite  in  line  with  the  expected  accuracies  for  this  study. 
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Reference  CSFs  were  determined  from  preliminary  Cl  calculations  at  180°  (linear)  and  140°.  The 
number  of  states  included  per  IRREP  are  3(At),  3(AU),  2(BU),  and  3(Bg),  and  the  contribution  of  the 
reference  set  to  the  wavcfunction  is  given  by  c)t!;s 

ref  CSFs 

r2  -  Y  r2 

i 

where  cf  is  the  CI  coefficient  for  the  i*  reference  CSF  in  a  given  state.  The  CI  results  at  both  angles 
give  cl„  £  (.88,  .90,  .90,  .88)  for  any  state  in  the  Ag,  Au,  Bu,  Bg  IRREPS,  respectively.  And,  for  the 
states  of  primary  interest,  cft/j  >  X  (.91),  2AU  (.91),  B  (.90).  The  list  of  reference  CSFs,  and  the  total 
number  of  CSFs  generated  in  each  IRREP,  are  shown  in  Table  2.  All  the  states  mentioned  above 
represent  electronic  excitations  out  of  the  In  Ky(bH)  In  ^(aj  orbitals.  The  excited  states  of  primary 
interest  are  obtained  by  excitations  into  the  In  ^(ap  lngl(bp  orbitals.  PESs  for  the  lowest  six  singlet 
states  are  shown  in  Figure  1,  and  the  data  given  explicitly  in  Table  3. 

3.1.2  Perturbed  Wavefunctions.  The  total  perturbed  wavefunctions,  'PfX)  and  4'(B),  arise  from 
the  magnetic  dipole  induced  perturbations  to  the  zeroth-order  electronic  state  wavefunctions,  VF°(X)  and 
^(B),  as  given  by 

¥(*)  =  ¥°(X)  +  'P'ffl/jX)  +  T‘(B/;X)  +  ^‘(A/iX)  03a) 

T(B)  -  ¥°(B)  +  'P1(A*;B)  +  ^‘(A^B)  +  ^¥l(B‘  J)  (13b) 


where  the  first-order  corrections  arise  from  the  magnetic  dipole  interactions, 

V'(Btx ;X)  :  <BJhJ'P0(X)> 
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Table  2.  Reference  CSFs  Used  in  Trans-bendng  Conformations  (C^  Symmetry)  and  Total  Number  of 
CSFs  Generated  Per  IRREP* 


MO 

IRREP 

2flg 

3as 

4a, 

5as 

2bu 

3bu 

4bu 

5bu 

6bu 

7bu 

la. 

lb. 

2bg 

Ag 

2 

2 

- 

- 

2 

2 

- 

- 

- 

- 

2 

- 

- 

(271428)b 

2 

2 

- 

2 

2 

- 

- 

- 

- 

- 

2 

- 

- 

2 

2 

- 

- 

2 

2 

- 

- 

- 

- 

- 

2 

- 

2 

2 

- 

1 

2 

1 

- 

- 

- 

- 

1 

1 

- 

2 

2 

- 

- 

2 

1 

1 

- 

- 

- 

2 

- 

- 

2 

2 

- 

- 

2 

1 

- 

- 

1 

- 

2 

- 

- 

2 

2 

- 

- 

2 

1 

- 

1 

- 

- 

2 

- 

- 

2 

2 

- 

- 

2 

1 

- 

- 

- 

1 

2 

- 

- 

K 

2 

2 

_ 

_ 

2 

1 

« 

_ 

• 

• 

2 

1 

(186914) 

2 

2 

1 

- 

2 

2 

- 

- 

- 

- 

1 

- 

- 

2 

2 

- 

1 

2 

2 

- 

- 

- 

- 

1 

- 

- 

2 

2 

- 

- 

2 

1 

- 

- 

- 

- 

2 

- 

1 

2 

2 

- 

- 

2 

1 

- 

- 

- 

- 

2 

1 

- 

2 

2 

1 

2 

2 

- 

- 

- 

- 

- 

1 

- 

- 

Bg 

2 

2 

_ 

2 

2 

1 

1 

_ 

_ 

(313910) 

2 

2 

- 

2 

2 

- 

1 

- 

- 

1 

- 

- 

2 

2 

- 

- 

2 

2 

- 

- 

1 

1 

- 

- 

2 

2 

- 

- 

2 

2 

- 

1 

- 

1 

- 

- 

2 

2 

- 

1 

2 

1 

1 

- 

- 

- 

1 

- 

2 

2 

- 

2 

2 

- 

- 

1 

- 

1 

- 

- 

2 

2 

- 

- 

2 

2 

- 

- 

- 

1 

1 

- 

- 

B, 

2 

2 

_ 

2 

2 

1 

1 

. 

(184782) 

2 

2 

- 

1 

2 

1 

- 

- 

- 

- 

2 

- 

- 

2 

2 

1 

- 

2 

1 

- 

- 

- 

- 

2 

- 

- 

2 

2 

- 

- 

2 

2 

- 

- 

- 

- 

1 

- 

1 

2 

2 

1 

- 

2 

1 

- 

- 

- 

- 

- 

2 

- 

■  The  Uj  and  lb,  MOs,  which  correspond  to  the  carbon  inner  shells,  were  treated  as  frozen  core  orbitals  in  the  Cl. 
b  Total  number  of  CSFs  actually  included  in  the  Cl  and  properties  calculations.  See  text  for  further  details. 
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Table  3.  Slate  Energies  vs.  Bend  Angle 


13 


I 


oo 
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:  <Bjpy  \^\X)> 
|p,  |¥0(X)> 
^GC*)  :  <i4j|i,|M»°(B)> 


V'^-B)  :  <A>y  |  ¥°(B> 


VHB^B)  :  <Bjp,  |¥°(B)> 

and  the  matrix  elements  listed  above  are  between  zeroth-order  wavefimctions  for  the  X  or  B  states  and 
the  manifold  of  states  having  the  specified  symmetry. 


3.1.3  %  Over  Perturbed  Wavefimctions.  Expanding  the  cross  products  and  dot  products  in 
Equation  3  and  inserting  the  first-order  corrections  to  the  wavefimctions  from  Equation  13  gives 

O 


Bf  = 


21m 

~T 


{ <'P°(X)  | my  1 4/0(B)  ><  X¥°(B)  | mt  |  T1  (B/;*)  > 


-  <  ^"(X)  | ro  |  ^(B)  ><  y¥°(B)  ImJ'P1  (B/;X)  > 

+  <'F°(X)  |  mz  |  *F°(B)><  T*(B')  |  my\  04/;X)> 

-  <  *P°(X)  \my  |4/0(B)><  4i0(B)  \mx  |  'P'(AJ,;X)>  } 


(14a) 


-  {<T0(X)  j  my  |  4,0(B)  ><  'F°(X>  |  m,  |  'P1(<;B)> 

-  <T°(X)  jmJ'f'0(B')x4'<>(X)  |m,  \'¥'(At;B)> 

*  <¥°(X)  |m,  !'P°(B')x'P,,(X)  \my  |'P1(B.‘;B)> 

-  <  ¥*(*)  |  my  |  y¥°(B)  ><  H^CX)  |  mx  \  ^'(B^J)  >}] 
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where  the  dagger  on  indicates  that  vibrational  motion  has  been  ignored  in  the  electronic 
0 

wavefunctions. 


Finally,  for  each  matrix  element  in  Equation  14.  we  need  to  obtain  the  vibrationally-averaged 
electric  transition  dipole  moment  between  two  bound  electronic  states.  For  example,  using  the  zeroth- 
order  PESs  to  calculate  the  vibrational  wavefunctions  we  get, 


S(X\B)lv:  =  <x?..(K)  \<^°(X)  \my  |'P°(fi)>  \%l ■(*)> 


(15a) 


5(B)  |  (B/;X) ):.,..  =  <xf  (R)  I  <V\B)  |  ffl,  |  'P1(B/;X)>  |  (*)>. 


(15b) 


The  vibrational  wavefunctions,  Xv  (R)’s,  are  obtained  by  numerically  solving  the  vibrational 
Schroedinger  equation  for  nuclear  motion  while  ignoring  rotational  effects.  The  integrals  in 
Equations  15a  and  15b  are  then  solved  numerically  using  Simpson’s  rule. 


The  final  form  for  which  includes  vibrational  averaging,  is 


B  = 


^L[{s(X\Byv.v..  •  5(B  |  (B/;X) )*.,..  -  S(X\B)'V.V..  •  S(B  |  (B/;X) ):., 
+  S(X  I  B)*,..  •  5(B  I  (A/;*)  )’ v..  -  S(X  |  B)J.,»  •  S(B_  |  (A/;X))*  v.) 

-  {5(X  | B-,.  •5(X|(A;;B'))'>..  -  5(X  |  B)*.y..  .5(^|(^;B));.V.. 

+  5(X  |  B)*,.  •  5(X  |  (B;;B')  -  5(X  |  ;£)?•.«  •  S(X  |  (Bj;B))*.v.j] 


(16) 


3.1.4  Quasi-degenerate  Perturbation  Theory.  The  B  and  2<\,  are  degenerate  at  180°  (and  nearly  so 
for  small  bend  angles);  thus,  the  magnetic  dipole  interaction  between  these  two  states  cannot  be  treated 
near  the  linear  conformation  using  the  non-degenerate  perturbation  theory  outlined  above.  Instead,  a 
quasi -degenerate  perturbation  scheme  (Yarkony  1988)  will  be  used  in  which  the  'FfAJ  space  is 
partitioned  into  two  parts;  a  P  space  containing  only  the  quasi-degenerate  electronic  state,  and 

a  space  containing  all  the  remaining  'F(A^)  states.  If  we  label  the  B  state  as  the  I’ th  state,  the 
wavefunction  for  B  becomes 
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¥'(/)  =  ¥«(/)  + 


(17) 


where  is  referred  to  as  a  dressed  wavefunction,  and  Tg  is  the  dressing  wavefunction.  'FgC/,/)  is 
defined  by 

^eC/,/)  *  GW,/);  withQ  «  (1  -  E  |  'P°(/)><  *¥°(J)  |)  (18) 

JeP(J) 

with  P(J)  containing  only  the  2\  element  for  the  current  case.  'FgC, /,/)  satisfies 

(H°  -  £)Tg(7./)  =  dpW)  (19) 

with  ii  the  magnetic  dipole  operator  and  Q  the  projection  operator  defined  in  Equation  18. 
Specifically,  &  projects  out  the  contribution  of  4/°(2AJ  from  'P'(J.I)  leaving  the  dressing  function, 
Tg (/,/),  which  contains  the  interaction  of  *F°(B)  with  all  other  *P°(AJ  except  the  'P°(2AU).  Equation 
19  is  then  solved  directly  in  the  CSF  basis  as  discussed  in  Section  2.2. 

We  now  want  to  include  the  interaction  of  the  quasi-degenerate  4/°(2Au)  with  the  *P°(B)  state  over 
the  magnetic  dipole  moment  operator,  p,  This  is  done  by  coupling  the  vibrational  levels  of  the  two 
electronic  states  to  form  new  vibronic  states,  'Flc(r,R)  (Yarkony  1988): 

Vtir.R)  -  XdZxlvfaR)  (20) 

IX 

where  I  runs  over  the  electronic  wavefunctions  for  the  two  quasi-degenerate  states.  The  Xx(£)  satisfy 

<f„  +  E,\R))xl(R)  -  xi(Ryt„  (21) 

with 

E,‘(R)  -  <'Vl(r-,R)\H'\'i’l(r,R)>. 
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The  *t'lc(r,R)  satisfy  the  total  (rotationless)  Schroedinger  equation, 


(Tn  +  H‘)^K(rJt)  *  EK^K(rJt), 


(22) 


where  K  is  a  coupled  vibronic  state  formed  as  a  linear  combination  of  the  vibrational  levels  of  the 
4>0(2AJ  and  VF°(B)  states  with  the  coupling  coefficients,  d,*,  determined  from  the  relationships 

(e  -  Ek)Dk  =  TDk  (23) 

with 


e/X./'X'  =  ^/,/^x,xe/x  ’ 

W-<x^li3 *'>*£>* 

-  <xilji(/./,)IXx->Jt  if  1*1' 

-  0  i/  /  *  /’ . 


(24) 


(25) 


The  magnitude  of  the  coupling  coefficients  is  determined  in  part  by  the  size  of  the  magnetic  dipole 
moment  matrix  element  coupling  the  2A^  and  B  states.  Based  on  the  coupling  coefficients,  dix,  one 
can  often  assign  the  vibronic  states,  (r,R),  as  being  composed  primarily  of  'F°(2AU)  or  VF°(B). 
Making  such  an  assignment  becomes  more  meaningful  when  |  eA  -  erx.  | »  T/x  rx. . 

The  contribution  to  the  3,-teim  in  Equation  16  through  the  matrix  elements  involving  the  vibronic 
wavefunctions,  4V(rJ0,  is  given  explicitly  by 

S(X  I  (a;,b))1.¥..  =  <xl (R)  \<^°(X)  I  m,  I  'v'q(a;-b)>  \ xi(R)> 

*  <  X*  (R)  I  <  V°(X)  |  m,  I  >  > .  (26) 
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From  Equation  20,  we  find  vFK(rJl)  to  be 


Vgirjt)  *  x  x?W'*'0W  *  XdfatWV'iB).  (27) 

X  Y 


Substituting  Equation  27  for  *FK  into  Equation  26  gives 

Stf  |  =  <%;.(/?)  |<^°(X)  |/«J'F1e(^;;fi)>  |xf(*» 

+  E4m*a<X?'(^)  |  <  'F'CX)  |mt  |’F°(2A11)>  |x?‘ “<*)> 

X  “  (28) 

♦  Y.dl,<X*-(R)  \<^°0C)  |ml'P°(B)>  \xUR)>- 

Y 


where  the  superscript  on  the  2\  in  d£»  indicates  the  coupling  coefficient  was  determined  using  the  x' 
component  of  the  magnetic  dipole  moment.  Based  on  symmetry  arguments,  the  third  term  in  Equation 
28  is  zero,  giving  the  final  expression 


S(X  |  -<  x?  (R)  |  <  ¥tf(X)  I  m,  |  V'eiAt  J)  >  \  Xl  (R)  > 

+  (R)  |<,F«(X)  |m,  |‘F°(2AJ>  |x?“(/0>. 

X 


(29) 


An  analogous  expression  can  be  written  for  the  S(X)  |  G4«;B)  )*.,..  matrix  element  in  Equation  16, 
with  the  coupling  coefficients,  A,  being  determined  from  the  y -component  of  the  magnetic  dipole 

moment  coupling  matrix  element 


3.2.1  Zeroth -order  State  Descriptions.  The  cis-bent  1B2  state  (in  Cjv  symmetry)  comes  from 
splitting  the  degenerate  1A„  and  is  described  (for  the  linear  molecule)  by  the  CSFs  (lo2  lo22o22o2) 

3o2  lji^ljt^lji^,  and  (.  .  ,)3o2IlJi2„lJt1„lJt1f,.  Its  degenerate  partner  is  the  linear  2A2  and  is  given 

by  the  CSFs  (.  .  .)3o2  l7t‘„l retire and  (.  .  OSo^lJt^lTt^lrc1,,.  The  active  orbitals  used  in  the 
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distributed  amongst  these  MOs.  The  2au(3b2)  MO  is  a  Rydberg  orbital.  Six  electronic  states  were 
then  averaged  including  lAlt  1,2B2,  1,2,3A2.  The  weighting  scheme  used  in  the  averaging  is 
wk=(l,  1,1, 1,1,1),  respectively.  The  SA-CASSCF  was  calculated  at  each  point  along  the  bending  PES, 
which  ranged  from  180°  (linear)  to  110°  at  10°  increments. 

Analogous  to  the  trans-bent  case,  the  reference  CSFs  were  determined  by  running  CIs  at  cis-bend 
angles  of  180°  (linear)  and  140°.  The  number  of  states  included  per  IRREP  were  3(A,),  3^,  2(B,), 
and  2(B2).  At  both  angles,  the  total  contribution  of  the  reference  set  to  the  wavefunctions  in  each 
IRREP  are  c)tfa  >  (.88,  .91,  .89,  .90),  respectively,  and,  for  the  states  of  primary  interest,  c*A  > 
X(A,)(.92),  1B2(.90).  The  list  of  reference  CSFs  and  the  total  number  of  CSFs  generated  in  each 
IRREP  are  shown  in  Table  4.  All  the  states  mentioned  above  represent  electronic  excitations  out  of 
the  l7Cux(b1)l7tuz(a1)  orbitals.  The  excited  state  of  primary  interest  here  is  the  1B2.  This  is  obtained 
by  excitations  into  the  lit^a^lTt^b*)  orbitals.  The  resulting  PESs  for  the  lowest  four  singlet  states 
are  shown  in  Figure  1  and  the  data  given  explicitly  in  Table  3. 

3.2.2  Perturbed  Wavefunctions.  The  perturbed  wavefunctions  'PCX)  and  'F(1B2)  may  be  written  as 

T(X)  «  VXX)  +  VXBfj)  +  T'(A2';X)  (30a) 

»F(1B2)  -  T-flBj)  +  ^(A*;  1 B2)  +  ¥*(£,*;  1 B2)  (30b> 

where  the  first-order  corrections  arise  from  the  magnetic  dipole  interactions 

'W;X)  :  <B2|pJ  ¥->(*)> 


W)  :  <A2lnI|'F»(X)> 

^\Ar,lB2)  :  <A,  |pJH'°(lB2)> 
ns,':  1*2)  :  <Bt  |p,  I'F°(1B2)> 

and  the  matrix  elements  listed  above  are  between  zeroth-ordcr  wavefunctions  for  the  X  or  1B2  states 
and  the  manifold  of  states  having  the  specified  symmetry. 


20 


Table  4.  Reference  CSFs  Used  in  Cl  for  Cis-bending  Conformations  C2v  Symmetry)  and  Total  Number  of  CSFs  Generated  Per  IRREP* 
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Tout  CSFs  per  IRREP. 


3.2.3  %  Over  Perturbed  Wavefunctions.  Expanding  the  cross  products  and  dot  products  in 
Equation  3  and  inserting  the  definitions  for  the  perturbed  wavefunctions  from  Equations  30  gives  ®0t 


b]  =  -^p[{<  ¥*tf)  I  m,  I  W  B2)  ><  B2)  j  mt  |  ^\B^X)  > 

-  <  r(X)  I  my  I  *P°(lfi2)><  ^(lB,)  I  mx  I  W/;*)} 

-  {  <  T°(X)  I  my  I  ><  ^\X)  I  m,  I  ^,'i  1B2)  > 

-  <^V\X)  \my  |4'°(1B2)x4'°(X)  |wJ'F,(B1,;l«a)>}].  (31) 

Finally,  for  each  matrix  element  in  Equation  31,  we  obtain  the  vibrationally-averaged  electric 
transition  dipole  moments  as  described  above. 

Using  the  nomenclature  established  for  the  B<— X  transition,  the  final  form  for  which  includes 
vibrational  averaging,  is 


Bo  -  I  ib2>;.,.-s(ib2  I  (B2X,X)K..-S(X  \  i B2y,¥..  • sob2  |  (A/;*)):.,.} 

-  {s(X  I  lB2yr.r. .  S(X)  I  (Af  ;  1B2»;.,.  - S(X  I  1B2),\.  • S(X  I  (B‘;  1B2))^}] .  (32) 


4.  RESULTS  AND  DISCUSSION 

4.1  B<-X  Transition;  Trans-bending.  Figure  2  shows  the  electric  dipole  transition  moments 
between  the  zeroth-order  states  including  <X\mM\B>,  <X\my\B>,  and  <X |mf  (24>i>  as  a  function 
of  bend  angle.  The  actual  values  are  listed  in  Table  5.  The  magnitudes  of  the  moments  are  essentially 
zero  at  180°  and  increase  with  bend  angle.  The  <X\my\B>  shows  a  maximum  near  140°,  which 
coincides  with  the  equilibrium  bend  angle  in  B  using  the  ground  state  bond  lengths.  Our 
<X|  mj  2A„>  curve  is  not  as  smoothly  varying  as  the  other  two  moments  in  Figure  2  due  to  an 
avoided  crossing  of  the  2A^  and  3A^  states  in  the  region  of  130-140°.  Also,  the  fact  that  <x|  mji> 
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<m>  laul 
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110  120  130  140  150  160  170 

Trans  Bend  Angle  (Degrees) 


180 

linear 


Figure  2.  Electric  Pinole  Transition  Moments  Between  the  Ground  State  and  the  Components  of  the 
l‘A..(2‘A...  1‘B-,  as  a  Function  of  Trans-bend  Angle.  All  Moments  are  Over  Zeroth-Order 
Wavefiinctions. 
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Table  5.  The  Electric  Dipole  Transition  Moments  (in  au)  vs.  the  Trans-bending  Angle  for  Matrix 
Elements  Involving  the  Zeroth-Order  Wavefiinctions  for  the  X  B,  and  2AU  States 


Angle  (Deg) 

<X\mx\B> 

<X\my\B> 

<X\mt\2A> 

180  (linear) 

-0.0394 

0. 

0. 

170 

0.0066 

160 

0.1071 

I 

150 

0.2224 

-0.4009 

-0.3283 

140 

0.3321 

-0.4276 

-0.4398 

130 

0.4294 

■ 

-0.6003 

120 

0.5189 

1 

-0.5790 

110 

0.5929 

-0.5686 

is  small  but  non-zero  at  180°  is  an  artifact  of  the  calculations  induced  by  a  minor  symmetry  breaking 
in  the  wavefunction  at  180°.  Its  effect  on  the  final  transition  intensities  was  tested  by  setting  the 
moment  to  zero  at  180°  and  rerunning  the  vibrational  analysis.  The  final  MCD  intensities  for  the  0-0 
and  0-1  transitions  changed  by  <  0.5%,  which  is  quite  acceptable  for  this  study. 


PPB  have  also  calculated  these  moments.  The  agreement  between  the  present  result  and  that  of 
PPB  is  quite  good  for  the  <x|  m,!  BE>  curve,  with  both  showing  a  maximum  near  140°.  Their 
<x|  mj  B>  and  <x|  m2 1 2A„>  moments  show  the  same  qualitative  dependence  on  bend  angle  as  what 
is  found  here.  However,  there  are  quantitative  discrepancies  between  the  studies  ranging  from  a 
few  percent  to  approximately  20%,  with  the  larger  difference  occurring  for  <x|mx|  B>  at  angles  less 
than  130°.  Given  the  differences  in  approach  between  the  two  studies,  including  their  use  of  a  smaller 
AO  basis  set,  different  bond  lengths,  and  a  different  Cl  approach,  it  would  be  very  difficult,  if  not 
impossible  (and  outside  the  purpose  of  this  study),  to  unambiguously  determine  the  source  of  these 
discrepancies  in  the  moments. 

4.1.1  Oscillator  Strengths.  The  absorption  oscillator  strengths  for  the  B<-X  transition  have  been 
calculated  and  are  given  in  Table  6.  The  1-0  band  carries  the  largest  oscillator  strength  of  ^=0.0025. 
This  can  be  compared  to  the  experimental  spectrum  of  Foo  and  Innes  (1973)  which  has  the  intensity 
peaking  around  v'=3-5.  It  is  difficult  to  make  any  further  comparison  with  experiment  due  to  the 
complicated,  diffuse  nature  of  the  spectrum,  and  the  absence  of  an  experimental  estimate  for  the 
oscillator  strength. 
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Table  6.  Calculated  Absorption  Oscillator  Strengths  for  the  B«-X  Transition  Resulting  From  the 
Vibrational  Analysis  on  the  Trans-bending  Mode 


v’-v" 

f^io3 

0-0 

0.905 

1-0 

2.490 

2-0 

0.478 

3-0 

0.281 

4-0 

0.301 

5-0 

0.146 

4.1.2  MCD  Intensities.  Equation  14  shows  the  transition  moments  needed  to  calculate  the  MCD 
for  the  trans-bending  mode.  There  are  two  main  parts  to  Equation  14.  The  first  part,  labeled  14a, 
accounts  for  the  perturbation  of  the  ground  state  and  the  second  part,  labeled  14b,  describes  the 
perturbation  of  the  excited  B  state.  Due  to  the  rather  large  energy  separation  of  the  ground  state  from 
the  nearest  excited  gerade  state  (i.e.,  in,:  Te£8.59eV  from  these  calculations),  one  would  expea  the 
first-order  corrections  to  X  to  be  considerably  smaller  than  the  corrections  to  B;  hence,  Equation  14a 
should  contribute  significantly  1ms  to  the  MCD  intensity  than  Equation  14b.  In  faa,  Segal  and  co- 
workers  (Dutch  and  Segal  1983;  Diamond  and  Segal  1984)  have  shown  this  to  be  true  in  other  cases, 
thus  often  enabling  one  to  ignore  the  first  term  in  Equation  14  at  a  considerable  computational  saving. 

Figure  3  shows  the  angle  dependence  of  the  orbital  angular  momentum  matrix  elements  between 
the  quasi-degenerate  states  B  and  2Au(<'F°(2Aj|Lx|'F°(it>  and  <'F°(2AU)  I  ly  I  T°(B)>)  over  their 
zeroth-order  wavefunctions.  For  a  ‘A  state  (linear  geometry),  the  matrix  element 
<4/°(2AJ  I L,  I  'F°(B)>  should  be  2  and  <M/0(2Ab)|L>|  T°(B)>  should  be  0,  consistent  with  the  results 
in  Figure  3.  The  <  *P°(2A11)  |  Lx  \  ^¥°{B)  >  value  drops  to  1.06  au  at  90°  while  the 
<  'F°(2 Au)  |  Ly  |  '¥\B)>  increases  to  0.75  au.  These  matrix  elements  are  necessary  for  determining 
the  perturbation  of  B  by  its  quasi-degenerate  partner  2A^.  Equations  16  and  29  show  the  contribution 
of  the  2AU-B  interaction  to  the  MCD  (with  v’=K)  to  be  proportional  to 

<  X?  (*)  I  <  ^tf)  K  I  >  I  *?•  <*)  ><  X*  GO  1  <  W0Z)  I  m,  |  VGA.) >  I E X^“ W  >  ■ 

X 
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Figure  3.  Orbital  Angular  Momentum  Matrix  Elements  <L>  Coupling  the  Zeroth-Order  Wavefunctions 
for  the  2'A.,  and  l‘Bu  Components  of  the  l‘Au.  The  Moments  are  Plotted  as  a  Function_of 
Trans-bend  Angle  (All  Values  are  in  Atomic  Units).  The  Magnetic^  Dipole  Moments  <m> 
are  Related  (in  au)  to  Orbital  Angular  Momentum  bv  <m>=('l/2')<L>. 
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The  PES  for  the  2A^  requires  some  discussion.  Peric,  Peyerimhoff,  and  Buenker  found  in  their 
calculations  that  the  2A„  state  prefers  a  non-planar  structure  with  a  twist  of  approximately  90°  about 
the  C-C  for  non-linear  conformations.  It  is  beyond  the  scope  of  this  study  to  attempt  a  multi¬ 
dimensional  vibrational  treatment,  therefore  the  molecule  is  assumed  to  remain  planar  in  the  2A„  state. 
This  limits  our  ability  to  calculate  an  accurate  vibronic  coupling  between  the  B  and  2\,  since  quasi¬ 
degenerate  perturbation  theory  requires  PESs  for  both  states.  In  analyzing  the  make-up  of  the  vibronic 
wavefunction,  VPK,  one  finds  the  encouraging  result  that  the  largest  contribution  to  the  first  two 
vibronic  states,  ¥,[^=0,1),  come  almost  entirely  from  couplings  with  the  v=0,l  of  the  2AU.  So  it  is 
reasonable  to  expect  the  *PK’s  to  be  good  approximations  to  the  true  wavefunctions  for  vK=0,l,  and  the 
calculations  to  predict  MCD  intensities  for  the  0-0  and  1-0  bands  in  reasonable  agreement  with 
experiment 

Figure  4  shows  the  plots  of  the  electric  dipole  transition  moment  matrix  elements  between  the 
zeroth-order  excited  B  state  and  the  first-order  correction  to  the  ground  state,  'F1(Bj;X).  These  matrix 
elements  are  found  in  Equation  14a.  Figure  5  plots  matrix  elements  between  the  zeroth-order  ground 
state  wavefunction  and  the  first-order  correction  to  the  B  state  (Equation  14b).  It  is  immediately  clear 
from  comparing  Figures  4  and  5  that  the  absolute  magnitudes  of  the  moments  involving  the  perturbed 
excited  B  state  are  considerably  larger  than  those  of  the  perturbed  X  state  over  most  of  the  PES. 

By  inspection  of  Figure  5,  one  sees  that  the  largest  matrix  element  in  the  vicinity  of  the 
equilibrium  bend  angle  for  the  B  state  is  <'F°(X)  |  mt  |  'PIC(A*;B)>,  which  represents  the  coupling  of 
B  with  the  A„  states  over  p*.  This  matrix  element  does  not  include  the  contribution  from  B  coupling  to 
its  quasi-degenerate  partner,  the  2A„.  It  is  left  to  the  vibrational  treatment  to  account  for  this  B-2AU 
interaction.  The  results  of  the  vibrational  treatment  are  given  in  Table  7,  which  lists  the  contributions 
from  each  matrix  element  in  Equation  14  to  each  vibrational  band.  The  matrix  elements 
<  '¥°(X)  |  mt  |  ¥’( A*\B) >  and  <  4/0(X)  |  mt  \  'V'iAfiB) >  now  contain  both  the  non-degenerate  and 
quasi-degenerate  (i.e.,  B-2AJ  PT  contributions.  It  can  be  seen  that  the  matrix  element  involving 
T'(A  *;B)  is  by  far  the  single  largest  contributor  to  the  intensity  for  the  two  vibrational  bands. 
Specifically,  the  B<— X  transition  gains  much  of  its  MCD  intensity  through  interactions  with  A„  states. 
As  in  absorption,  the  more  intense  MCD  band  is  the  1-0.  The  calculated  are  (0-0)=-3.48, 

(l-0)=-5.82.  Table  7  also  lists  (in  parenthesis)  the  corresponding  experimental  Ae““values; 

(0-0)=-0.8,  (1-0)=- 1.6.  The  theoretical  values  for  (0-0)  and  (1-0)  are  a  factor  of  4  larger  than  the 
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Figure  4.  Electric  Dipole  Transition  Moments  as  a  Function  of  Trans-bend  Angle.  The  Matrix 

Elements  Include  the  First-Order  Correction  to  the  Wavefunction  for  the  Perturbed  Ground 
State.  All  Values  are  in  Atomic  Units.  See  Equation  13  for  Definitions  of  First- Order 
Wavefunctions. 
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Figure  5.  Electric  Dipole  Transition  Moments  as  a  Function  of  Trans-bend  Angle  for  &— X  The 
Matrix  Elements  Include  the  Total  First-Order  Correction  to  the  Wavefunction  for  the 
Perturbed  B  State.  All  Values  are  in  Atomic  Units.  See  Equation  13  for  Definitions  of 
First-Order  Wavefunctions. 
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Table  7.  Electric  J)ipole  Transition  Moment  Components  (in  au)  Determining  the  MCD  Intensity  for 
the  B<— X  Transition  (Trans-bending  Motion) 


0-0 

v'-v" 

1-0 

<VP°(X)  \mx  |'P°(B)> 

-.0201149 

+.0105815 

Zeroth-order 

<¥•(*)  |mj  ¥#(B)> 

+.0696649 

-.1190800 

Moments 

<'P°(B)|mJ'P1(fi;;X)> 

-.0180523 

+.0280062 

<'F°(B)|mJ'P1(B/;X)> 

+.0175951 

-.0354671 

Perturbing  tne 

<r(B)|mjl'Fl(A;;i)> 

+.0114457 

-.0566385 

X  state 

<'F°(B)|mJ'F1(4;;X)> 

-.0246575 

+0.3196930 

<4'°(X)|  mt\^\AxuJ)> 

+16.3245000 

-19.7214400 

<'F°(X)|mJ'P1(^;B)> 

-.8072640 

-2.7946300 

Perturbing  the 

<4/0(X)  \my  \ 

-1.6048100 

-4.9070150 

B  state 

<^°(X)\mx  m*:;*)> 

+.1523540 

+.2918810 

V 

-.224X104 

-.463x10-* 

-3.48 

-5.82 

MCD  Band  Maxima 

Ae™"b(cgs  units) 

(-0.8)c 

(-1.6) 

(Experimental  Maxima) 

*  Consistent  with  electric  dipole  moments  in  debye. 

b  Gaussian  band  shape  described  in  text  with  A=120cm  '(0-0),  150ciri'(l-0),  and  14,  =54116cm  ')0-0),  54794cm  ‘(l-0). 
e  Experimental  values  from  Gedenken  and  Schnepp  (1976). 


experimental  maxima,  but  of  the  correct  sign.  This  agreement  is  acceptable  considering  the 
approximate  nature  of  the  vibrational  treatment. 

4.2  1B,+-X  Transition:  Cis-bending.  Theory  suggests  (vide  supra)  that  the  1B2<-X  transition 
might  overlap  the  B+-X  band  system  and  be  responsible  for  part  of  the  irregular  band  shapes  and 
variations  in  intensities.  In  this  section,  we  examine  the  predicted  MCD  and  absorption  intensities 
(i.e.,  oscillator  strengths)  for  the  1B2<-X  transition  and  compare  these  to  the  analogous  values  for  the 
B«-X  system. 
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4.2.1  Oscillator  Strengths.  The  y-component  of  the  electric  dipole  transition  moment  for  1B2<— X 
is  shown  as  a  function  of  bending  in  Figure  6.  Again,  the  1B2  state  is  one  of  the  degenerate  partners 
in  the  1 state  (see  Figure  1);  thus,  the  electric  dipole  transition  moment  must  be  zero  at  180°.  The 
curve  in  Figure  6  shows  the  moment  continually  increasing  from  essentially  zero  at  180°  to 
0.58  au  at  110°,  the  largest  angle  calculated.  Once  again,  Peric,  Peyerimhoff,  and  Buenker  have 
calculated  this  electric  dipole  moment  and  found  a  similar  qualitative  dependence  on  bend  angle. 
However,  the  magnitude  of  thechange  per  degree  is  larger  in  PPBs  study  than  here.  For  example,  at 
an  angle  of  1 10°  they  predict  a  moment  of  =  1 0.80 1  while  our  value  is  1 0.58 1 .  The  reader  is  referred 
to  the  section  concerning  the  B<— X  transition  for  a  discussion  of  possible  causes  for  this  quantitative 
difference  in  the  moments  between  the  two  studies.  Table  8  lists  the  values  of  the  electric  dipole 
transition  moment  plotted  in  Figure  6. 

The  calculated  0-0  transition  energies  for  the  B<— X  and  1B2<— X  band  systems  are 
AE=56667cm'1(B«-X)  and  AE=56330cm'1(lB2<— X),  which  do  indeed  lie  close  to  one  another.  The 
calculated  transition  energies  for  excitations  from  v"=0  to  v'=0,l,2,3  are  given  in  Table  9  along  with 
the  experimental  values  for  the  B«-X  reported  by  Foo  and  Innes  (1973). 

Table  8.  Electric  Dipole  Transition  Moment  vs.  Cis-bend  Angle  for  the  Matrix  Element  Involving  the 
Zeroth-order  X  and  1B2  States 


Angle  (Deg) 

<X(1'A,) \my  1 11B2> 

180  (linear) 

■ 

170 

160 

0.1461 

150 

0.2662 

140 

130 

0.4659 

120 

0.5327 

110 

0.5798 
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Table  9.  Transition  Energies  (in  cm"1)  for  the  B*-X  and  1B2<— X  Band  Systems 


v'-v" 

B<-X* 

A  E(s'm-v\) 

B<-5? 

A  E(yM-v\) 

1B2<-X 

.  1 

56667. 

54116. 

56330. 

57296. 

629. 

54850. 

734. 

56988. 

57756. 

460. 

55538. 

688. 

57585. 

1 

58398. 

642. 

56285. 

743. 

58070. 

4-0 

59227. 

829. 

56984. 

699. 

58523. 

5-0 

60170. 

943. 

57827. 

843. 

59116. 

*  From  PESs  using  minimum  energy  band  lengths  for  both  states.  See  text.  Section  3. 
b  Experimental  values  from  Foo  and  Innes  (1973). 


The  oscillator  strengths  for  absorption  are  given  in  Table  10.  The  transitions  to  the  first  four 
v'  levels  are  shown,  and  one  finds  the  2-0  band  predicted  to  be  the  most  intense  with  ^=4.4x10 s,  and 
15  times  more  intense  that  the  0-0.  The  important  point  to  note  is  that  the  magnitude  of  the  2-0  is 
approximately  a  factor  of  60  less  intense  than  the  strongest  band  in  the  Bt-X 


Table  10.  Calculated  Absorption  Oscillator  Strengths  for  the  1B2<-X  Transition 


v'-v" 

f^io5 

04) 

0.2820 

1-0 

1.6800 

2-0 

4.3800 

3-0 

3.0200 

4-0 

0.0146 

4.2.2  MCD  Intensities.  Unlike  the  B  state,  the  1B2  has  no  allowed  magnetic  dipole  coupling  to  its 
quasi-degenerate  partner,  the  2A^.  This  is  significant  because  the  theoretical  study  of  Peric, 
Peyerimhoff,  and  Buenker  found  the  2Aj  flike  the  2AJ  component  of  the  1A„  prefers  a  twisted 
conformation  at  larger  bend  angles.  The  lack  of  magnetic  dipole  coupling  between  1B2  and  2A2 
simplifies  the  calculation  of  the  MCD  for  the  1B2<— X  because  quasi -degenerate  PT  is  not  needed.  But, 
more  importantly,  this  removes  the  uncertainty  in  the  calculated  MCD  added  by  having  to  couple  to  a 
PES  (i.e.,  the  2Aj)  which  is  known  to  be  inaccurate  at  higher  vibrational  levels. 
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Table  1 1  gives  the  MCD  Ae““  for  the  first  six  vibrational  bands.  The  vibrational  levels  4-0  and 
5-0  are  found  to  be  opposite  in  sign  to  the  first  four  bands.  One  important  observation  from  Table  1 1 
is  that  matrix  elements  involving  the  perturbed  ground  state  contribute  significantly  to  the  total  MCD 
intensity  for  the  0-0,  1-0,  and  2-0  bands.  This  arises  not  because  of  large  perturbation  to  the  ground 
state,  but  because  the  two  components  describing  the  perturbation  to  the  1B2  state  are  nearly  equal  in 
magnitude  and  cancel  one  another.  This  is  in  contrast  to  what  we  saw  in  the  B«-X  band  system  where 
essentially  all  the  MCD  intensity  arose  from  perturbations  of  the  excited  state  B,  not  the  ground  state. 
The  other  important  observation  is  that  the  strongest  band,  the  3-0,  has  a  predicted  MCD  intensity 
Ae”“=  -.0037 ,  which  is  three  orders  of  magnitude  smaller  than  in  the  B<— X  Part  of  the  exP^anat'on 
for  this  relatively  small  MCD  signal  lies  in  the  symmetry  of  the  molecule  which  dictates  the  allowed 
magnetic  dipole  couplings  and  therefore  the  ability  to  borrow  transition  intensity.  It  will  be 
informative  to  look  at  the  states  closest  in  energy  to  the  lA^  since  these  should  couple  most  strongly 
with  the  components  of  the  1A„. 

Peric,  Buenker,  and  Peyerimhoff  (1984,  1987)  calculated  the  lowest  16  singlet  states  in  the  linear 
conformation,  and  the  energies  of  the  lowest  five  states  are  shown  in  Table  12,  with  the  spread  in  AE 
due  to  different  choices  in  reference  CSFs.  As  a  reminder,  the  1A„  splits  into  the  (2A2,  lB^)  (cis- 
bending)  and  (2A,,  1BU)  (trans-bending).  Table  12  shows  its  nearest  neighbors  to  be  the  ligand  lfl^. 

It  has  already  been  demonstrated  that  in  the  trans-bending  case,  the  B  state  can  borrow  intensity 
from  any  A„  and  B„  state.  This  includes  its  quasi-degenerate  partner,  the  2A^,  and  the  components 
making  up  both  the  1Z^  (AJ  and  lll^A^,  BJ  states.  Since  these  four  neighboring  "states"  (in  lie 
within  *0.35e  of  the  B  state  (in  the  linear  geometry),  perturbation  theory  would  predict  potentially 
large  borrowing  of  intensity  from  these  four  states  by  B,  and  hence  a  strong  MCD  signal. 

If  the  same  rationale  is  applied  to  the  cis-bending  mode,  its  Qv  symmetry  disallows  any  interaction 
of  the  1B2  state  with  any  Aj  state,  including  its  quasi -degenerate  partner  the  2A2  and  its  lower  energy 
neighbor  the  1Z^  (A^.  Therefore,  in  the  reduced  Cjv  symmetry,  the  1B2  can  borrow  intensity  from 
only  two  of  its  four  closest  neighbors  which  are  the  components  of  the  llTJA,,  B,).  The  next  nearest 
allowed  coupling  is  with  the  1A,,  which  lies  «1.3eV  (linear  geometry)  above  the  1A„.  Peric,  Buenker, 
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Table  1 1.  Electric  Dipole  Transition  Moment  Components  (in  au)  Determining  the  MCD  Intensity  for  the  1B2<— X  Transition 
(Cis-bending  Motion)* 
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Table  12.  Veitical  Excitation  Energies  for  Linear  Acetylene* 


State  (DoJ 

AE(eV) 

in. 

8.38-8.46 

in. 

7.92-8.00 

ia. 

7.56-7.65 

i£; 

7.21-7.32 

is; 

0 

*  From  ab  initio  MRD-CI  calculations  (Peric,  Buenker,  and 
Peyerimhoff  1984). 


and  Peyerimhoff  (1984,  1987)  also  reported  the  cis-  and  trans-PESs  associated  with  the  16  lowest 
singlet  D,^  states,  and  these  curves  clearly  show  the  AE  between  the  1B2  and  all  higher  lying  states 
increases  rapidly  with  bend  angle.  In  fact,  their  curves  predict  that  this  increase  is  more  pronounced 
for  the  1B2  than  for  B,  again  suggesting  less  ability  for  the  1B2  to  borrow  intensity  from  neighboring 
states  than  what  we  see  for  the  B  state. 

Therefore,  both  perturbative  (i.e.,  based  on  energetics)  and  symmetry  arguments  can  be  used  to 
explain  the  calculated  results  which  show  the  MCD  intensity  to  be  significantly  weaker  in  the  1B2<— X 
transition  than  in  the  Bt-X 

5.  CONCLUSIONS 

A  new  method  has  been  presented  for  predicting  a  priori  the  MCD  spectrum  in  molecules.  The 
method,  based  on  first-order  perturbation  theory,  significantly  reduces  the  time  consuming  process  of 
explicitly  calculating  a  large  number  of  excited  electronic  states  which  might  be  needed  to  satisfy  the 
summation  in  the  first-order  perturbation  expression  in  Equation  9.  As  an  initial  test  case,  this  method 
has  been  used  here  to  study  the  character  of  the  experimental  MCD  spectrum  associated  with  the  Bt-X 
electronic  transition  in  acetylene.  The  MCD  and  absorption  intensities  for  the  1B2<— X  transition  were 
also  determined  to  assess  the  possibility  that  the  irregularities  seen  in  the  185-170  run  spectral  region 
result  from  the  1B2«— Xband  system  overlapping  the  Bt-X 
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For  the  B*-X  band  system,  this  study  predicts  a  Ae"“  of  -3.5  and  -5.8  for  the  0-0  and  1-0 
vibrational  bands,  respectively.  These  results  agree  with  experiment  with  respect  to  the  signs  of  the 
bands.  They  are  also  consistent  with  the  experimentally  observed  trend  of  the  Ae”“  increasing  as  v' 
increases.  However,  the  magnitudes  are  a  factor  of  4  larger  than  the  experimental  values  of  -0.8 
and  -1.6,  respectively.  This  agreement  is  quite  acceptable  considering  the  approximations  made  in  the 
vibrational  treatment. 

Concerning  the  possible  overlap  of  the  1B2«— X  and  B<— X  band  systems,  the  current  study  predicts 
the  0-0  transition  energies  for  the  two  systems  to  lie  in  the  same  spectral  region.  The  calculated  0-0 
transition  energies  are  AE=56330cm'1(lB2<— X)  and  AE-56667cm1(B«— X).  The  most  intense  MCD 
band  in  the  1B2<— X  band  system  is  the  3-0  with  a  calculated  A -.0037.  For  the  B«— X  band 
system,  the  more  intense  vibrational  band  is  predicted  to  have  a  value  of  A£“(l-0)*  -5.82,  which  is 
-103  times  larger  than  the  Ae^n  values  calculated  for  the  1B2«— X  MCD  spectrum.  In  spite  of  the 
possibility  that  improvements  in  the  theoretical  treatment,  especially  the  vibrational  treatment,  might 
cause  significant  changes  in  the  calculated  intensities,  it  seems  unlikely  that  the  MCD  intensities  for 
the  1B2<-X  and  B«-X  would  change  three  orders  of  magnitude  relative  to  one  another.  Hence,  this 
study  predicts  that  the  irregularities  in  the  185-170  nm  region  of  the  MCD  are  not  due  to  the 
overlapping  of  the  1B2<— X  and  B<— X  band  systems. 

With  regards  to  the  absorption  intensity,  the  current  study  predicts  that  the  oscillator  strengths  for 
the  lBjt-X  transition  are  a  factor  of  60  smaller  than  those  calculated  for  the  B<— X  suggesting  that  the 
1B24-X  would  at  most  appear  as  small  perturbations  in  the  absorption  spectrum  of  B<— X.  One 
interference  in  the  Bt-X  absorption  spectra  has  already  been  established  by  Foo  and  Innes  (1973)  who 

M  mm 

identified  high-lying,  trans-bending  vibrational  levels  of  the  AOAJt-X  transition  spilling  over  into  the 
B*-X,  with  the  A*-X(13-0)  occurring  at  an  energy  slightly  above  the  0-0  for  B<— X.  Another  possible 
candidate  is  the  2AV«-X,  which  is  found  here  and  by  previous  theory  (Peric,  Buenker,  and  Peyerimhoff 
1987)  to  be  strongly  electric  dipole  allowed  at  larger  bend  angles  and  a  planar  conformation. 

However,  at  these  bend  angles,  the  molecule  was  found  to  prefer  a  twisted  conformation,  leading  one 
to  expect  a  poor  Franck-Condon  overlap  with  the  ground  state;  and  thus,  predicting  a  potentially  weak 
transition  (Peric,  Buenker,  and  Peyerimhoff  1987).  And  finally,  as  Peric,  Peyerimhoff,  and  Buenker 
had  mentioned,  there  are  most  likely  strong  vibrational  couplings  of  the  2B„  state  with  other  states.  In 
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light  of  the  results  from  this  study  and  previous  work  by  others,  these  couplings,  along  with  the  At-X 
band  system,  seem  to  be  the  most  plausible  explanations  for  the  diffuse  bands  and  complicated 
structure  in  the  absorption  spectrum  of  the  Ek— X  band  system. 

In  summary,  this  study  predicts  that  the  1B2«-X  and  B*-X  electronic  transitions  should  occur  in 
the  same  spectral  region.  However,  these  calculations  do  not  support  the  hypothesis  that  the  lB^X 
band  system  makes  a  significant  contribution  to  either  the  absorption  or  the  MCD  spectrum  in  the 
185-170  nm  region. 
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1  Commandant 

USAFAS 

ATTN:  ATSF-TSM-CN 
Fort  Sill,  OK  73503-5600 

1  F.J.  Seiler 

ATTN:  S.  A.  Shackleford 
USAF  Academy,  CO  80840-6528 

1  University  of  Dayton  Research 

Institute 

ATTN:  D.  Campbell 
AL/PAP 

Edwards  AFB,  CA  93523 
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1  NASA 

Langley  Research  Center 
Langley  Station 

ATTN:  G.  B.  Northam/MS  168 
Hampton,  VA  23365 

4  National  Bureau  of  Standards 

ATTN:  J.  Hastie 
M.  Jacox 
T.  Kashiwagi 
H.  Semerjian 

US  Department  of  Commerce 
Washington,  DC  20234 

1  Aerojet  Solid  Propulsion  Co. 

ATTN:  P.  Micheli 
Sacramento,  GA  95813 

1  Applied  Combustion  Technology, 

Inc. 

ATTN:  A.  M.  Varney 
P.O.  Box  607885 
Orlando,  FL  32860 

2  Applied  Mechanics  Reviews 

The  American  Society  of  Mechnical 
Engineers 
ATTN:  R.  E.  White 
A.  B.  Wenzel 
345  E.  47th  Street 
New  York,  NY  10017 

1  Atlantic  Research  Corp. 

ATTN:  M.  K.  King 
5390  Cherokee  Avenue 
Alexandria,  VA  22314 

1  Atlantic  Research  Corp. 

ATTN:  R.  H.  W.  Waesche 
7511  Wellington  Road 
Gainesville,  VA  22065 

1  AVCO  Everett  Research 

Laboratory  Division 
ATTN:  D.  Stickler 
2385  Revere  Beach  Parkway 
Everett,  MA  02149 

1  Battelle  Memorial  Institute 

Tactical  Technolgoy  Center 
ATTN:  J.  Huggins 
505  King  Avenue 
Columbus,  OH  43201 


No.  of 

Copies  Organization 

1  Cohen  Professional  Services 
ATTN:  N.  S.  Cohen 
141  Channing  Street 
Redlands,  CA  92373 

1  Exxon  Research  &  Eng.  Co. 

ATTN:  A.  Dean 
Route  22E 

Annandale,  NJ  08801 

1  Ford  Aerospace  and 

Communications  Corp. 
DIVAD  Division 
Div.  Hq.,  Irvine 
ATTN:  D.  Williams 
Main  Street  &  Ford  Road 
Newport  Beach,  CA  92663 

1  General  Applied  Science 

Laboratories,  Inc. 

77  Raynor  Avenue 
Ronkonkama,  NY  11779-6649 

1  General  Electric  Ordnance 

Systems 

ATTN:  J.  Mandzy 
100  Plastics  Avenue 
Pittsfield,  MA  01203 

2  General  Motors  Rsch  Labs 
Physics  Department 
ATTN:  T.  Sloan 

p  Tootc 

Warren,  Ml  48090 

2  Hercules,  Inc. 

Allegheny  Ballistics  Lab. 
ATm:  W.  B.  Walkup 
E.  A.  Yount 
P.O.  Box  210 

Rocket  Center,  WV  26726 

1  Alliant  Techsystems,  Inc. 

Marine  Systems  Group 
ATTN:  D.  E.  Braden/ 

MS  MN50-2000 
600  2nd  Street  NE 
Hopkins,  MN  55343 

1  Alliant  Techsystems,  Inc. 

ATTN:  R.  E.  Tompkins 
MN38-3300 
5700  Smetana  Drive 
Minnetonka,  MN  55343 
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1  IBM  Corporation 

ATTN:  A.  C.  Tam 
Research  Division 
5600  Cottle  Road 
San  Jose,  CA  95193 

1  IIT  Research  Institute 

ATTN:  R.  F.  Remaly 
10  West  35th  Street 
Chicago,  IL  60616 

2  Director 

Lawrence  Livermore  National 
Laboratory 

ATTN:  C.  Westbrook 

M.  Costantino 
P.O.  Box  808 
Livermore,  CA  94550 

1  Lockheed  Missiles  &  Space  Co. 

ATTN:  George  Lo 
3251  Hanover  Street 
Dept.  52-35/B204/2 
Palo  Alto,  CA  94304 

1  Los  Alamos  National  Laboratory 

ATTN:  B.  Nicholas 

T 7,  MS-B284 
P.O.  Box  1663 
Los  Alamos,  NM  87545 

1  National  Science  Foundation 

ATTN:  A.  B.  Harvey 
Washington,  DC  20550 

1  Olin  Ordnance 

ATTN:  V.  McDonald,  Library 
P.O.  Box  222 

St.  Marks,  FL  32355-0222 

1  Paul  Gough  Associates,  Inc. 
ATTN:  P.  S.  Gough 

1048  South  Street 
Portsmouth,  NH  03801-5423 

2  Princeton  Combustion 

Research  Laboratories,  Inc. 
ATTN:  M.  Summerfield 

N.  A.  Messina 
475  US  Highway  One 
Monmouth  Junction,  NJ  08852 

1  Hughes  Aircraft  Company 

ATTN:  T.  E.  Ward 
8433  Faibrook  Ward 
Canoga  Park,  CA  91303 


No.  of 

Copies  Organization 

1  Rockwell  International  Corp. 
Rocketdyne  Division 
ATTN:  J.  E.  Flanagan/HB02 
6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

4  Sandia  National  Laboratories 

Division  8354 
ATTN:  R.  Cattolica 
S.  Johnston 
P.  Mattern 
D.  Stephenson 
Livermore,  CA  94550 

1  Science  Applications,  Inc. 

ATTN:  R.  B.  Edelman 
23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

3  SRI  International 

ATTN:  G.  Smith 
D.  Crosley 
D.  Golden 

333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  Stevens  Institute  of  Tech. 

Davidson  Laboratory 
ATTN:  R.  McAlevy,  III 
Hoboken,  NJ  07030 

1  Sverdrup  Technology,  Inc. 

LERC  Group 

ATTN:  R.  X  Locke,  MS  SVR-2 
2001  Aerospace  Parkway 
Brook  Park,  OH  44142 

1  Thiokol  Corporation 

Elkton  Division 
ATTN:  S.  F.  Palopoli 
P.O.  Box  241 
Elkton,  MD  21921 

1  Thiokol  Corporation 

Huntsville  Division 
ATTN:  J.  Deur 
Huntsville,  AL  35807-7501 

3  Thiokol  Corporation 

Wasatch  Division 
ATTN:  S.  J.  Bennett 
P.O.  Box  524 
Brigham  City,  UT  84302 
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United  Technologies  Research 
Center 

ATTN:  A.  C.  Eckbreth 
East  Hartford,  CT  06108 

United  Technologies  Corp. 
Chemical  Systems  Division 
ATTN:  R.  S.  Brown 

T.  D.  Myers  (2  copies) 
P.O.  Box  49028 
San  Jose,  CA  95161-9028 

Universal  Propulsion  Company 
ATTN:  H.  J.  McSpadden 
Black  Canyon  Stage  1 
Box  1140 

Phoenix,  AZ  85029 

Veritay  Technology,  Inc. 

ATTN:  E.  B.  Fisher 
4845  Millersport  Highway 
P.O.  Box  305 

East  Amherst,  NY  14051-0305 

Brigham  Young  University 
Dept,  of  Chemical  Engineering 
ATTN:  M.  W.  Beckstead 
Provo,  UT  84058 

California  Institute  of  Tech. 

Jet  Propulsion  Laboratory 
ATTN:  L.  Strand,  MS  512/102 
4800  Oak  Grove  Drive 
Pasadena  CA,  91109 

California  Institute  of  Technology 
ATTN:  F.  E.  C.  Culick/MC  301-46 
204  Karman  Lada. 

Pasadena,  CA  91125 

University  of  California 
Los  Alamos  Scientific  Lab. 

P.O.  Box  1663,  Mail  Stop  B216 
Los  Alamos,  NM  87545 

University  of  California,  Berkeley 
Chemistry  Department 
ATTN:  C.  Bradley  Moore 
211  Lewis  Hall 
Berkeley,  CA  94720 

University  of  Calfornia,  San  Diego 
ATTN:  F.  A.  Williams 
AMES,  B010 
La  Jolla,  CA  92093 


2  University  of  California, 

Santa  Barbara 
Quantum  Institute 
ATTN:  K.  Schofield 
M.  Steinberg 

Santa  Barbara,  CA  93106 

1  University  of  Colorado  at  Boulder 
Engineering  Center 

ATTN:  J.  Daily 
Campus  Box  427 
Boulder,  CO  80309-0427 

2  University  of  Southern  California 

Dept,  of  Chemistry 

ATTN:  S.  Benson 
C.  Wittig 

Los  Angeles,  CA  90007 

1  Cornell  University 

Department  of  Chemistry 
ATTN:  T.  A.  Cool 
Baker  Laboratory 
Ithaca,  NY  14853 

1  University  of  Delaware 

ATTN:  T.  Brill 
Chemistry  Department 
Newark,  DE  19711 

1  University  of  Florida 

Dept,  of  Chemistry 
ATTN:  J.  Winefordner 
Gainesville,  FL  32611 

3  Georgia  Institute  of  Technology 
School  of  Aerospace  Engineering 
ATTN:  E.  Price 

W.  C.  Strahle 
B.  T.  2nn 
Atlanta,  GA  30332 

1  University  of  Illinois 

Dept,  of  Mech.  Eng. 

ATTN:  H.  Krier 

144MEB,  1206  W.  Green  Street 
Urbana,  IL  61801 


No.  of 


No.  of 

Copies  Organization 


Johns  Hopkins  University/APL 
Chemical  Propulsion  Information 
Agency 

ATTN:  T.  W.  Christian 
Johns  Hopkins  Road 
Laurel,  MD  20707 

University  of  Michigan 
Gas  Dynamics  Lab 
Aerospace  Engineering  Lab 
ATTN:  G.  M.  Faeth 
Ann  Arbor,  Ml  48109-2140 

University  of  Minnesota 
Dept,  of  Mechanical  Engineering 
ATTN:  E.  Fletcher 
Minneapolis,  MN  55455 

Pennsylvania  State  University 
Applied  Research  Laboratory 
ATTN:  K.  K.  Kuo 

H.  Palmer 
M.  Micci 

University  Park,  PA  16802 

Pennsylvania  State  University 
Dept,  of  Mechanical  Engineering 
ATTN:  V.  Yang 
University  Park,  PA  16802 

Polytechnic  Institute  of  NY 
Graduate  Center 
ATTN:  S.  Lederman 
Route  110 

Farmingdale,  NY  11735 

Princeton  University 
Forrestal  Campus  Library 
ATTN:  K.  Brezinsky 

I.  Glassman 
P.O.  Box  710 
Princton,  NJ  08540 


1  Purdue  University 

School  of  Aeronautics  and 
Astronautics 
ATTN:  J.  R.  Osborn 
Grissom  Hall 

West  Lafayette,  IN  47906 


1  Purdue  University 

Department  of  Chemistry 
ATtN:  E.  Grant 
West  Lafayette,  IN  47906 


2  Purdue  University 

School  of  Mechanical  Engineering 
ATTN:  N.  M.  Laurendeau 
S.  N.  B.  Murthy 
TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 

1  Rensselaer  Polytechnic  Inst. 

Dept,  of  Chemical  Engineering 
ATTN:  A.  Fontijn 
Troy,  NY  12181 

1  Stanford  University 

Dept,  of  Mechanical  Engineering 
ATTN:  R.  Hanson 
Stanford,  CA  94305 

1  University  of  Texas 

Dept,  of  Chemistry 
ATtN:  W.  Gardiner 
Austin,  TX  78712 

1  University  of  Utah 

Dept,  of  Chemical  Engineering 

ATTN:  G.  Flandro 

Salt  Lake  City,  UT  84112 

1  Virginia  Polytechnic  Institute  and 

State  University 
ATTN:  J.  A.  Schetz 
Blacksburg,  VA  24061 

1  Freedman  Associates 

ATTN:  E.  Freedman 
2411  Diana  Road 
Baltimore,  MD  21209-1525 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3204 _ Date  of  Report  February  1991 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


r 

5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 

,  .  •* 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
Changes  to  organization,  technical  content,  format,  etc.)  . . 


"  •*?.*•* 

Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground.  MD  2101-5066 

OFFICIAL  BUSINESS 


NO  POSTAGE 
NECESSARY 
IF  HALED 
M  THE 

UNTTED  STATES 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001,  APG,  MD 


POSTAGE  WIU.  BE  PAID  BY  ADDRESSEE 


Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T- 

Aberdeen  Proving  Ground,  MD  21005-9989 


-FOLD  HERE 


